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We study the properties of chaos in the motions of a spinning test particle in 
Schwarzschild spacetime. We characterize the chaos using the power spectrum of the 
time series of z components of the particle's position. It is found that the pattern 
of the power spectrum shows not only white noise but also l//-type fluctuation, 
depending on the value of the total angular momentum J and the spin S of the 
test particle. Therefore we succeed in classifying the chaotic motions, which have 
been classified as simply chaotic ones in former works, into the two distinct types. 
One is 1//, and the other is white noise. Based on this classification, we plot, in 
the two-dimensional parameter space (J, S), the phase diagram for the properties of 
the chaos. This phase diagram enables us in principle to guess the properties of the 
system (J, S) by observing the dynamics of the test particle, even if the motion is 
chaotic. Furthermore, we detect that the origin of the 1/f fluctuation is that the 
particle motion stagnates around regular orbits (tori), while traveling back and forth 
between them, which is called "stagnant motion" or "sticky motion" in Hamiltonian 
dynamical systems. The point is that the difference of the property of the chaos or 
the power spectra is due to the topological structure of the phase space, which in 
turn is governed by the physical parameter set (J, S) of the system. From this point 
of view, the chaos we found in this system is not always merely random. 
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I. INTRODUCTION 



Nature is filled with phenomena that exhibit chaotic behavior. In chaotic systems we 



cannot predict the system's future state exact 
found in some relativistic systems js, Q, 0, [? 



Such chaotic behavior has a. 
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2l|. For example, in Schwarzschild spacetime, the motions of a spinning test particle can be 



chaotic 



13| . If the test particle does not have spin, the motion of the test particle is regular. 



They have found that, as the magnitude of the spin increases extremely, the motions switch 
from regular to chaotic, using the Poincare map and the Lyapunov exponent. In practice, 
the magnitude of the spin where the particle motion is remarkably chaotic is so large that 
such a system is not realistic, which has been remarked in the paper [l^. However, this 
model is important in understanding chaos in general relativistic systems. 

In this paper, we look for statistical laws to characterize the chaos in the motions of 
the spinning test particle in Schwarzschild spacetime. Actually, the chaotic motions in this 
system have been classified merely as chaotic, according to the distribution of the points in 
Poincare maps and positiveness of the Lyapunov exponents, but the details of the properties 
of the chaos have not been clarified 131] • Once we find chaotic behavior, however, we should 
rather characterize it to extract the specific properties of the system. Indeed, we can hardly 
learn anything about the chaos if we judge it only from the randomness of the distribution 
of the points in Poincare maps or the positiveness of the Lyapunov exponents l|, [2I, . Not a 
few people believe that a chaotic system is simply random and completely unpredictable. It 
is true that we cannot predict the time evolution of the state of the system exactly, when the 
system is chaotic. However, we should note that, even in such cases, we can frequently find 
some statistical laws which are specific to the system. One possible measure of chaos is the 
power spectrum of the time series of the system. If the power spectrum is white noise, the 
time evolution of the system is not time-correlated. In not a few cases, however, the pattern 



of the power spectrum obeys power laws, so-called 1// fluctuations 22 
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251,1261,1271, 



which can be clearly distinguished from the white-noise type. That the power spectrum 
obeys some power laws means that the time evolution of the system is time-correlated. 
Thus we can classify the chaos from the pattern of its power spectrum. 

Let us explain our strategy to characterize the chaos in the motions of the spinning test 
particle in Schwarzschild spacetime in this paper. To begin with, we introduce the power 
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spectrum of the time series of z components of the particle's position. Next we characterize 
the properties of the chaos, using the pattern of the power spectrum. It is found that 
the pattern of the power spectrum can be classified as 1// or white noise. That is, we 
succeed in classifying the chaotic motions into two distinct types. Furthermore, based on 
this classification, we can plot, in the two-dimensional parameter space {J,S), the phase 
diagram for the properties of the chaos. Finally, we detect the origin of the 1// fluctuations 
of the power spectrum in this system. We find out that the orbit stagnates around the 
tori, while traveling back and forth between them, whenever the power spectrum shows 1// 
spectral pattern. 

This paper is organized as follows. In Sec. [Ill we shall briefly review the basic equations, 
i.e., the equations of motion for a spinning test particle in Schwarzschild spacetime. In 
Sec. mil we plot the Poincare maps of the chaotic motions in this system. Here we point out 
a weakness in the Poincare maps, which brings a motivation to introduce another method 
to characterize such chaotic motions. In Sec. IIVI we introduce the power spectrum to char- 
acterize the properties of the chaos in this system. Then we find the pattern of the power 
spectrum can be classified into two types, 1// and white noise, depending on the values 
of the system parameters. In Sec. |V] we detect the origin of the 1// fluctuations. The 
final section is devoted to summary and discussion. Throughout this paper we use units 
c = G=l 



II. EQUATIONS FOR A SPINNING TEST PARTICLE IN SCHWARZSCHILD 

SPACETIME 



We consider a spinning test particle in Schwarzschild spacetime. 



ds' = - 1 



2M 



2M 



dr^ + r'^dO'^ + sin' 



(1) 



where M is the mass of the black hole. The equations of motions of a spinning test particle 
in relativistic spacetime have been derived by Papapetrou 28|] and then reformulated by 



Dixon 



29| . The set of equations is given as 

dx^ 



dr 

Dpf" 1 



dr 



(2) 
(3) 
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= 2pK''\ (4) 

dr 

where r, v^,p^, and S'^'^ are an affine parameter of the orbit, the four- velocity of a particle, 
the momentum, and the spin tensor, respectively, deviates from a geodesic due to the 
coupling of the Riemann tensor with the spin tensor. We adopt the additional condition 



formulated by Dixon 



29|, 



P.S^'' = 0, (5) 
which gives a relation between and v^, and consistently determines the center of mass of 
the spinning particle. The mass of the particle /i is defined by 

/i' = -PmP^ (6) 

To make clear the freedom of this system, we have to check the conserved quantities. Re- 
gardless of the symmetry of the background spacetime, it is easy to show that the mass fi 
and the magnitude of spin S defined by 

^2 = ^S^^S'^" (7) 

are constants of motion [sojl . If a background spacetime possesses some symmetry described 
by a Killing vector 



is also conserved 



Q = - l^,;uS'''' (8) 

2^. Because the spacetime we consider in this paper is static and spher- 



E = = -Pt - (9) 



ically symmetric, there are two Killing vector fields, and ^'^^y From ([8]), we find the 
constants of motion related with those Killing vectors as 

M 

J, = C(<^) =p^- r(5'^" - r5^^ cot 6) sin^ 6. (10) 

E and are interpreted as the energy of the particle and the z component of the total 
angular momentum, respectively. Because the spacetime is spherically symmetric, the x 
and y components of the total angular momentum are also conserved. In addition, without 
loss of generality we can choose the z axis in the direction of total angular momentum as 

(Jx.,J„J.) = (0,0,J), (11) 

where J > 0. In the following sections, we integrate the above equations of motion nu- 
merically for various values of parameters E, J, and S, using the Bulirsch-Stoer methods 



3l|. 



5 



III. POINCARE MAPS 



In this section, we illustrate the particle motions in the model formulated in the previous 



section (see also 
the paper 
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28 



291]) using Poincare maps as shown in Fig. [H As mentioned in 
13| . the parameter range of the energy E for each fixed parameter set (J, S) is 
restricted to enable the particle to move around the black hole without going to infinity or 
'ailing into the black hole. Although similar analyses have been already done in the paper 
isl . the analysis in this section motivates us to introduce in the next section another method 
to characterize the motions. 

Using Poincare map, we can judge if the motions are chaotic or not. To plot the Poincare 
map, we adopt the equatorial plane {6 = vr/2) as a Poincare section and plot the point (r, 
v^) when the particle crosses the Poincare section with < (Fig. [1]). In Fig. [21 we plot 
the Poincare maps for the total angular momentum J = AfiM. The values of the spin S 
are set to S = 1.2/iM in Fig. [2] (a) and S = 1.4/iM in Fig. [2] (b). Each value of the energy 
E is chosen as an appropriate one so that the test particle does not escape to infinity and 
does not fall into a black hole. The parameter sets (J, S, E) in Figs. [2] (a) and [2] (b) are 
the same as those in Figs. 4 (e) and 4 (f) in the paper y. In Fig. [2] dots with different 
colors correspond to the data from the orbits with different initial conditions. If the orbit 
is chaotic, some of the tori are broken and the Poincare map no longer consists of a set 
of closed curves. Both in Fig. [2] (a) and Fig. [2] (b), the points in the Poincare maps are 
scattered randomly, and the so-called chaotic sea is formed. That is, the orbits with both 
sets of parameters are chaotic. 

Now it is worthwhile to note a weakness in the Poincare map. The chaotic sea in Figs. [2] 
(a) and[2](b) cannot be distinguished apparently. Indeed, the existence of chaotic sea in the 
Poincare map allows us to judge whether the particle motions in this model are regular or 
chaotic {13I]. The positiveness of the Lyapunov exponent, which for this model have been 
also investigated in jisj], does too. That is, we can tell that the motion is chaotic using these 
measures. However, only from these measures, we cannot know more than that the motion 
is merely chaotic, when it is chaotic. Once we find the chaotic behaviors, we should rather 
characterize the chaos, since being chaotic does not always mean randomness or no rule. 
Therefore, another method is necessary to characterize the chaos in more detail. In the next 
section, we will introduce the power spectrum to classify such chaotic motions. 
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IV. 1// FLUCTUATIONS OF THE POWER SPECTRUM 

In this section, we characterize the chaos in the spinning test particle motions which was 
shown in the previous section. Here we analyze the time series of the particle position. In 
order to do that, first of all, we introduce the power spectrum. The power spectrum of the 
time series of z components of the particle's position, Pziy^)-, is defined by 



P. feu) 







(12) 



Here we set T = 10^ in our computation. If we define the autocorrelation function ^^(t) as 

$.(r) = ^ I ^z{t)z{t + T)dt, (13) 



and take the limit T oo, Wiener-Khinchin's theorem 
Py(uj) and the autocorrelation function ^,(t) as 



32j relates the power spectrum 



oo 



P^icu) = / <I>,(r)e-^-^rfr. (14) 



oo 



We plot the power spectrum Pziyj) in Fig. [31 In Figs. [3] (a) and[3](b), we choose parameter 
sets with the same values as those in Figs. [2] (a) and [2] (b), respectively. Each line color 
in Figs. [3] (a) and [3] (b) also corresponds to that of the dots in Figs. [2] (a) and [2] (b), 
respectively. Note that the patterns of the power spectra are different between Fig. |3] (a) 
and Fig. [3](b). In particular, we find that one is 1//*^ (/ = oj/2'k) where z/ ~ 1.2 (Fig. [3](a)), 
while the other is white noise (Fig. [3](b)) in the low-frequency range about uj < 0.01. Using 
Eq. f|T^ . we see that if the power spectrum is of l/f type, the temporal correlation is also 
of the power-law type, which means a strong temporal correlation with no typical time scale. 
On the other hand, if the power spectrum is of white noise type, temporal correlation is a 
5-function, which means no temporal correlation. Therefore we can clearly distinguish the 
chaotic motions by using the pattern of the power spectrum. Moreover, Fig [3] shows that 
the pattern of the power spectrum Pz{uj), that is the property of the chaos, depends on the 
system parameters, the total angular momentum J and the spin S of the test particle. In 
addition, as shown in Fig [3], the pattern of the power spectrum is independent of the initial 
conditions, if we fix the parameter set. 

Furthermore, we test the pattern of the power spectrum Pz{uj) for various grid points in 
the two-dimensional space (J, S), and obtain a new phase diagram for the properties of the 
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chaos, summarized in Fig. HI Each value of the energy E is chosen as an appropriate one so 
that the test particle does not escape to infinity and does not fall into the black hole, and 
the chaotic motions that we pay attention to in this paper are dominant. In the blank region 
in Fig. m the test particle goes to infinity or falls into a black hole, since the energy surface 
is unbounded. This blank region agrees with the region of the type (U2) in Fig. 3 in the 
paper In the region where circle (Q) cind cross (x) symbols are marked in Fig. IH the 
particle motion is chaotic and bounded. In the region where circle symbols (Q) are marked 
in Fig. m the power spectrum Pz{oj) shows l//-type spectrum. On the other hand, in the 
region where cross symbols (x) are marked, the power spectrum Pz{uj) shows white noise. 
The total region together with dotted and crossed regions is included in the region of the 
type (B2) in Fig.3 in the paper p^]. In the paper jl3], the motions in the region of the type 
(B2) are classified as to be merely chaotic ones. Therefore, Fig. H] means that we succeed in 
classifying the chaotic motions for various values of parameter set (J, S), which have been 
classified as simply chaotic in former works, into two distinct types, 1/f and white noise, 
using the power spectrum Pz{uj). 

Now we consider the physical meanings of the tendency of the power spectrum pattern 
to change from 1// to white noise by increasing the magnitude of spin S of the test particle 
in Fig. m Let us note that the system is integrable and the motion of the test particle is 

nn 

not chaotic, if the test particle does not have spin (5 = 0) [2|,113[. Furthermore, th 
of the test particle with the small magnitude of spin remains almost regular 33 
The motion become remarkably chaotic at last, if the magnitude of the spin is extremely 
large [3] • Then it is reasonable to regard the magnitude of the spin S of the test particle as 
a measure of chaos. Considering these facts, our results can be interpreted to mean that the 
pattern of the power spectrum Pz{uj) changes from 1/f type to white noise as the strength 
of the chaos increases. Therefore, Fig. H] means that we succeed in classifying the chaotic 
motions into two categories in accordance with the strength of chaos. In the next section, 
we investigate the origin of the 1/f fiuctuations in the power spectrum Pz{uj). 

V. ORIGIN OF THE 1/f FLUCTUATIONS IN THE POWER SPECTRUM 



34 



35| 



In this section, we investigate the origin of the 1/f fiuctuations in the power spectrum 
Pz{uj) that we found out in the previous section. The relations between the physical quanti- 
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ties and figures are summarized in Table HI To begin with, we look more closely at the time 
evolution of the orbits which were analyzed in the previous section (see Fig. [3]). 

First, in Fig. [5l we plot the time series of the components in the Poincare map with 
z = and < 0. We choose parameter sets in Figs. [5] (a) and [5] (b) with the same values 
as those in Figs. [2] (a) and [2] (b), respectively. In Fig. [5] (a) the values of the components 
in the Poincare map fluctuate with stagnation. It seems that the duration of the stagnation 
varies from time to time, and is a mixture of various lengths. There is no typical one. This 
feature is observed commonly when the power spectrum Pz{uj) is the 1// type. Then we 
expect that this behavior is closely related to the 1/f fluctuations in the power spectrum 
Pz{u)) in Fig. [3] (a). In Fig. [5] (b), on the other hand, the values of the v"^ components in 
the Poincare map oscillate almost monotonously. This behavior is consistent with the power 
spectrum Pz{uj) showing white noise in Fig. [3t^b). 

Second, in Fig. [6l we plot the same orbit we shown in Fig. [5] (a) in the two-dimensional 
configuration space {r,z). We find that the whole orbit (Fig. [6] (a)) has two significant 
components (Figs. M (b) and [6] (c)). The periods when the orbit stagnates around each 
component in Fig. [6] perfectly correspond to that when components in the Poincare map 
stagnate in Fig. [5] (a). We never see such stagnant motions when the power spectrum 
become white noise (Fig. EJ^b)). These results suggest that the 1/f fiuctuations in the 
power spectrum Pz{uj) are originated by the stagnation and the stickiness around these two 
distinguishable components, while traveling back and forth between them, in the particle 
motions. 

Now, let us show a mechanism by which the motion is stagnant or not, as shown in 
Figs. [5] and [6l that is, the power spectrum Pz{uj) shows 1/f or white noise as shown in 
Figs [3|, depending on the combination of {J,S), as shown in Fig. HI The essence of the 
mechanism is that such stagnant motions are originated around the tori (regular orbits), 
and the structure of the tori changes depending on the combination of {J,S). To see the 
structure of tori more clearly, it is useful to see the motions at lower energy where the chaotic 
sea shrinks and tori get larger. In Fig. [7| we plot the Poincare maps with same values of J 
and 5* as in Fig. [2] but smaller values of E. Around this value of energy E, the structure of 
tori is not very sensitive to the change of E, compared to the changes of J or S. The volume 
of the chaotic sea in the Poincare maps in Fig. [7] is smaller than that in Fig. [2l However, we 
will see that the topological structure of the phase space is quite similar to that of Fig. O 
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First, we pay attention to the chaotic orbits which we saw in Fig. [71 In Fig. [8] we plot 
the time series of the v"^ components in the Poincare maps Fig. [3 We adopt each data 
in Figs. M (a) and M (b) as the chaotic orbits that are plotted in Figs. [7| (a) and [7] (b), 
respectively. It is found that the values of the components fluctuate with stagnation in 
Fig. M (a), while they do not in Fig. M (b). The feature about the stagnation in Figs. [7] (a) 
and[7] (b) is similar to that in Figs. [5] (a) and[5] (b), respectively. In Fig. [9l we plot the same 
orbit we showed in Fig. [8] (a) in the two-dimensional configuration space {r,z). Similarly 
to the higher-energy case in Figs. [5] and [6l the periods when the orbit stagnates around 
each component in Fig. [9] perfectly correspond to that when components in the Poincare 
map stagnate in Fig. [8](a). Now, in Fig. [10] we plot the power spectrum Pz{uj) of the same 
chaotic orbits as Figs. M (a) and [HI (b). It is clearly confirmed that the pattern of the power 
spectrum Pz{uj) shows 1// in Fig. [10] (a), while the pattern shows white noise in Fig. [TO] (b). 
Therefore, similarly to the higher-energy case, it is suggested that the 1/f fluctuations in 
the power spectrum Pz{oj) are originated by the motion which stagnates and sticks around 
the two distinguishable components (Figs. [9] (b) and [9] (c)), while traveling back and forth 
between them. Our results indicate that the pattern of the power spectrum Pz{uj) of the 
chaotic orbits strongly depends on the value of the total angular momentum J and spin S, 
while the value of energy E does not affect the pattern of the power spectrum so strongly, 
although it affects the volume of the tori. 

Second, we pay attention to the tori (regular orbits) which we saw in Fig. [7] In Fig. [TT] 
we plot the orbits of the tori in the two-dimensional configuration space (r, z) . Each color 
of lines in Figs. [TT] (a) and [H] (b) corresponds to that of closed curves in Figs. [7] (a) and [7] 
(b). For example, the orbits A and B in Fig. [TT](a) correspond to the closed curves A and B 
in Fig. [7] (a). Indeed, the existence of some regular orbits in this model have been shown in 
the paper [l3|. In this paper, however, we should emphasize that the structure of the phase 
space is quite different between panels (a) and (b) in Fig. [7] or Fig. [TTJ For example, any 
tori which correspond to A and B in Fig. [7] (a) or Fig. [TT] (a) can never be found in Fig. [7] 

(b) or Fig. [TT] (b). That is, we should emphasize that the topological structure of the tori 
changes depending on the parameter sets (J, S), and the change is in accordance with the 
change of types of the power spectra (see Figs. [TO] and [TT]) . Moreover, it is quite important 
to note that the orbits of the tori A and B in Fig. [TT] (a) resemble the chaotic orbits in Fig. [6] 

(c) and Fig. [6] (b), or Fig. [9] (b) and Fig. [9] (c), respectively. This resemblance means that 
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the orbits in Figs. [6] and [9] stagnates around the tori A and B in Fig. [TT] (a), although it is 
not easy to detect such tori in Fig. [2] (a) since chaotic sea dominates and the tori A and B 
shrink extremely with the larger value of the energy E. 

Summing up the above results leads us to the following conclusions. The 1/f fluctuations 
in the power spectrum are originated by the orbit which stagnates around the tori while 
traveling back and forth between them. Whether we can observe 1/f fluctuations in the 
power spectrum in the time series of a component depends on the topological structure of 
the phase space. It is the pair of tori A and B in Figs [7] (a) and [11] (a) that gives rise to 
the 1/f fluctuations in the power spectrum Pz{uj). Moreover, our results suggest that the 
topological structure of the tori changes strongly depending on the total angular momentum 
J and spin S rather than the energy E. Therefore, it is reasonable that the pattern of power 
spectrum Pz{oj) of the chaotic motion depends mainly on the parameter sets (J, S), and we 
can expect that the phase diagram that we obtained in the previous section (see Fig. H]) is 
reliable. 



VI. SUMMARY AND DISCUSSION 



In this paper we have characterized the properties of chaos in a spinning test particle 
in Schwarzschild spacetime. We have calculated the power spectrum of the time series of 
z components of the test particle's position, Pz{uj), and found out that the pattern of the 
power spectrum Pz{uj) is 1/f or white noise in the low- frequency range (see Fig. [3]). That 
is, we have succeeded in classifying the chaotic motions, which had been classified as merely 
chaotic in the paper [l^ , into these two distinct types (see Table ^ . The important point 
is that the pattern of the power spectrum strongly depends on the spin S and the total 
angular momentum J of the test particle and not on the initial conditions. Our analyses 
also suggest that the value of the energy E does not affect so strongly whether the power 
spectrum Pz{uj) becomes l//-type spectrum or not (see Figs. [3] and [T0|) . Then, testing the 
pattern of Pz{uj) for various grid points in two-dimensional {J, S) plane, we have obtained 
the phase diagram for the character of the chaotic motions (see Fig. Hj). This phase diagram 
enables us in principle to guess the properties of the system (J and S) by observing the 
dynamics of the test particle, even if the motion is chaotic. 

Furthermore we have pointed out that the pair of tori A and B in Figs [7] (a) and [H] 
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(a) gives rise to the l//-type power spectrum of the time series of z components of the 
particle's position PziiS). The important point is that the chaotic orbits stick near the tori, 
while traveling back and forth between them. Whenever the power spectrum Pziio) becomes 
the 1// spectral pattern, we have found that the orbit (Fig. [6] (a)) stagnates around two 
significant components (Figs. [6] (b) and [6] (c)). Moreover, investigating the motions with 
the value of lower energy where more tori dominate, we have confirmed that such significant 
components imitate the orbits of the tori characterized as A and B in Figs. [7] (a) andfTT] (a). 

Eventually, the conclusion is summarized as follows. We have two types of chaos as we 
have seen by the power spectra in the system we have studied in this paper. One is 1//, 
and the other is white noise. The difference of the properties of chaos or the power spectra 
is caused by the topological structure of the phase space, which in turn is governed by the 
physical parameter set (J, S*) of the system. From this point of view, the chaos we found in 
this paper is not always merely random. 

The type of motion where the phase point in chaotic orbit stays close to some regular 
orbits (tori) for some long time is known as "stagnant motion" or "sticky motion," and is 



often observed in Hamiltonian dynamical systems 36|, l37|, l38|| . Stagnant motions are usually 



accompanied by 1/ f fluctuations and are considered to be due to the fractal structure of 



the phase space 39|, |40|, |41|. In particular, stagnant motions are often observed for weakly 



chaotic, nearly integrable systems. This is consistent with the fact that we have observed 
the 1// fluctuations for smaller values of spin S when the chaos is weak, and white noise 
for larger values of spin S when the chaos is strong. Until now such 1// fluctuations have 
not been discovered in any relativistic systems. We have shown that the 1// fluctuations 
we observed for the first time in the relativistic system are also generated by such stagnant 
motion, that the particle motion stagnates around regular orbits, while traveling back and 
forth between them in Schwarzschild spacetime. 

Finally the astrophysical implications of our results should be mentioned. It is true that 
the value of the spin S where the particle motion is remarkably chaotic is so large that such 
a star cannot exist, which has been remarked in the paper [isl. Then the result obtained in 
this paper is relevant as an illustration and is not directly applicable to a physical system. 
However, it is expected that the method and theory in this paper can be useful and applied 
to other astrophysical systems in practice, since nature is filled with phenomena that exhibit 
chaotic behavior. 



12 



Acknowledgments 

We would like to thank Kei-ichi Maeda for valuable discussions. H.K. and K.K. are sup- 
ported by JSPS Fellowship for Young Scientists. The authors are grateful to the anonymous 
referee for helpful comments on improving the text and the figures. 



[1] A. J. Lichtenberg and M. A. Lieberman, "Regular and Chaotic Dynamics", (Springer 1992). 

[2] G. Contopoulos, "Order and Chaos in Dynamical Astronomy", (Springer 2002). 

[3] Deterministic Chaos in General Relativity, edited by D. Hobill, A. Burd, and A. Coley 

(Plenum, New York, 1994), and references therein. 

[4] J. D. Barrow, Phys. Rep. 85, 1 (1982). 

[5] G. Contopoulos, Proc. R. Soc. London A431, 183 (1990). 

[6] C. P. Dettmann, N. E. Prankel and N. J. Cornish, Phys. Rev. D50, R618 (1994). 

[7] U. Yurtsever, Phys. Rev. D52, 3176 (1995). 

[8] V. Karas and D. Vokrouhhcky, Gen. Rcla. Grav. 24, 729 (1992). 

[9] H. Varvoglis and D. Papadopoulos, Astron. Astrophys. 261, 664 (1992). 

[10] L. BombeUi and E. Calrzetta, Class. Quantum Grav. 9, 2573 (1992). 

[11] R. Moeckel, Commun. Math. Phys. 150, 415 (1992). 

[12] Y. Sota, S. Suzuki, and K. Maeda, Class. Quantum Grav. 13, 1241 (1996). 

[13] S. Suzuki and K. Maeda, Phys. Rev. D55, 4848 (1997). 

[14] S. Suzuki and K. Maeda, Phys. Rev. D58, 023005 (1998). 

[15] S. Suzuki and K. Maeda, Phys. Rev. D61, 024005 (1999). 

[16] M. D. Hartl, Phys. Rev D67, 024005 (2003). 

[17] O. Semerak, Mon. Not. R. Astron. Soc. 308, 863 (1999). 

[18] J. Levin, Phys. Rev. Lett. 84, 3515 (2000). 

[19] J. D. Schnittman and F. A. Rasio, Phys. Rev. Lett. 87, 121101 (2001). 

[20] N. J. Cornish and J. Levin, Phys. Rev. Lett. 89, 179001 (2002). 

[21] K. Kiuchi and K. Maeda, Phys. Rev. D70, 064036 (2004). 

[22] P. Dutta and P. M. Horn, Rev. Mod. Phys. 53, 497 (1981). 

[23] M. B. Weissman, Rev. Mod. Phys. 60, 537 (1988). 



13 



[24] H. G. Schuster, "Deterministic Chaos" (Physik-Verlag, 1984) 
[25] P. Manneville, Journal de Physique, 41, 1235 (1980). 
[26] T. Kohyama, Prog. Theor. Phys. 71, 1104 (1984). 

[27] T. Geisel, A. Zacherl and G. Radons, Phys. Rev. Lett. 59, 2503 (1987). 
[28] A. Papapetrou, Proc. Roy Soc. A209, 248 (1951). 

[29] W. G. Dixon, Proc. Roy. Soc. A314, 499 (1970); ibid A319, 509 (1970); Journ. Gen. Rel. 
Grav. 4, 193 (1973); Phil. Trans. Roy. Soc. 277, 59 (1974); in Isolated Gravitating Systems in 
General Relativity, edited by J. Ehlers, (North-Holland, Amsterdam, 1979), 156. 

[30] E. Wald, Phys. Rev. D 6, 406 (1972). 

[31] W. Press, B. P. Flannery, S. Teukolsky, and W. T. Vetterling, Numerical Recipies in C (Cam- 
bridge University Press, Cambridge, England, 1986). 
[32] L. E. Leichl, A Modern Course in Statistical Physics (Univ. of Texas press, 1980) 
[33] A. N. Kolmogolov, Dokl. Akad. Nauk. SSSR 98, 527 (1954). 
[34] V. I. Arnold, Russ. Math. Surveys 18, 85 (1963). 
[35] J. Moser, Nachr. Akad. Wiss. Gottingen. Math. Phys.,1, 1 (1962). 
[36] C. F. F. Karney, Physica D8, 360 (1983). 

[37] B. V. Chirikov and D. L. Shepelyansky, Physica D13, 395 (1984). 

[38] Y. Aizawa, Y. Kikuehi, T. Harayama, K. Yamamoto, M. Ota and K. Tanaka, Prog. Theor. 

Phys. Suppl. 98, 36 (1989). 
[39] Y. Aizawa, Prog. Theor. Phys. 71, 1419 (1984). 
[40] J. D. Meiss and E. Ott, Phys. Rev. Lett. 55, 2741 (1985). 
[41] Y. Y. Yamaguchi and T. Konishi, Prog, of Theor. Phys. 99, 139 (1998). 



14 



TABLE I: summary of figures 



s 


1.2nM 


1.4;uM 




l/r (Figs. El (a), do] (a)) 


white noise (Figs. [3] (b), M (b)) 




stagnating (Figs. [5] (a), [8] (a)) 


non stagnating (Figs. [5] (b), [8] (b)) 


Tori A, B 


observed (Figs. El [7] (a) , [3 E] (a) ) 


not observed (Figs. [3 (b), E] (b)) 



TABLE II: Scliematic classification of motions in this system 



bounded 



regular • • • (periodic or quasi periodic) 



chaotic 



1//^ • • • (correlated) 

white noise ■ • • (uncorrelated) 



unbounded 
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FIG. 1: The schematic drawing of the Poincare maps. 

FIG. 2: The Poincare maps with z = and < 0. All orbits have the total angular momentum 
J = 4.0/uM. The magnitude of spin and the total energy are S = 1.2fj,M and E = 0.93545565^ in 
panel (a), and S = 1.4/iM and E = 0.92292941// in panel (b), respectively. The dots with different 
colors correspond to the data from the orbits with different initial conditions. 

FIG. 4: The phase diagram for the type of the power spectrum pattern. The patterns of the power 
spectrum for the chaotic orbits at grid points in a two-dimensional (J, S) configuration are tested. 
At the points where circle symbols (Q) ai^e marked, the l//-type power spectrum is observed. At 
the points where cross symbols (x) are marked, the white-noise power spectrum is observed. At 
the points where triangular symbols (A) are marked, the orbit apparently behaves almost regular 
in the temporal interval of numerical computation. In the blank region, the test particle goes to 
infinity or falls into the black hole, since the energy surface is unbounded. 

FIG. 5: The time series of the v''"{t) components of the Poincare map with z = and < 0. 
Parameter sets (J, S, E) in panels (a) and (b) are chosen so as to have the same values as those in 
Figs. [2] (a) and (b), respectively. 

FIG. 6: The orbit which we saw in Fig. (a) in the two-dimensional configuration space (r, z). (a) 
The orbit for the whole period < t < 50000. (b) The orbit for the period 4000 < t < 11000. (c) 
The orbit for the period 24000 < t < 28000. 

FIG. 7: The Poincare maps with z = and < 0. Parameter sets (J, S) in panels (a) and (b) are 
chosen as the same values as those in Figs. [2] (a) and [2] (b), respectively. The values of the total 
energy are chosen as E = 0.933/x (panel (a)) and E = 0.9205/i (panel (b)), that are smaller than 
that of Fig. [5] (a) and [5] (b), respectively. 

FIG. 3: The power spectrum of the time series of z components of the particle. Each set of the 
parameters and the initial conditions in panels (a) and (b) are the same as in Figs. [2] (a) and 
(b), respectively. The long-time correlations with the power law, so called 1/f fluctuations, are 
observed in panel (a), while no such correlations are observed in panel (b). 
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FIG. 8: The time series of the v'''[t) components of the Poincare map with z = and v < . Each 
orbit corresponds to that of the chaotic orbits in Figs. [7] (a) and[7](b), respectively. 



FIG. 9: The orbit which we saw in Fig. [8] (a) in the two-dimensional configuration space (r, z). (a) 
The orbit for the whole period < t < 70000. (b) The orbit for the period 22000 < t < 38000. (c) 
The orbit for the period 42000 < t < 70000. 



FIG. 10: The power spectrum Pz{uj) of the chaotic orbits in Figs [7] (a) and[7](b). We adopt the 
time-series z{t) in panels (a) and (b) as the chaotic orbits in Figs. [7] (a) and [7] (b), respectively. 
The long-time correlations with power law, so called 1/f fluctuations, are observed in panel (a), 
while no such correlations are observed in panel (b). 



FIG. 11: The orbits in the two-dimensional configuration space (r, z) corresponding to each torus 
in Figs. [7] (a) and [7] (b). Each line color in panels (a) and (b) corresponds to that of the closed 
curves in Figs. [7] (a) and[7](b), respectively. For example, the regular orbits A and B in Fig. [9] (a) 
correspond to the tori A and B in Fig. [7] (a). 
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